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ABSTRACT 

It has long been accepted that a possible mechanism for explaining the existence 
of magnetic white dwarfs is the merger of a binary white dwarf system, as there 
are viable mechanisms for producing sustainable magnetic fields within the merger 
product. However, the lack of rapid rotators in the magnetic white dwarf population 
has been always considered a problematic issue of this scenario. Smoothed Particle 
Hydrodynamics simulations show that in mergers in which the two white dwarfs have 
different masses a disc around the central compact object is formed. If the central 
object is magnetized it can interact with the disc through its magnctosphere. The 
torque applied by the disc changes the spin of the star, whereas the transferred angular 
momentum from the star to the disc determines the properties of the disc. In this work 
we build a model for the disc evolution under the effect of magnetic accretion, and 
for the angular momentum evolution of the star, which can be compared with the 
observations. Our model predicts that the magnetospheric interaction of magnetic 
white dwarfs with their discs results in a significant spin down, and we show that for 
magnetic white dwarfs with relatively strong fields (larger than 10 MG) the observed 
rotation periods of the magnetic white dwarf population can be reproduced. We also 
investigate whether turbulence can be sustained during the late phases of the evolution 
of the system. When a critical temperature below which turbulence is not sustained is 
introduced into the model, the periods of the three fast rotating, strongly magnetic, 
massive white dwarfs in the solar neighborhood are recovered. 

Key words: stars: white dwarfs, stars: magnetic field, accretion discs 



1 INTRODUCTION of massive white dwarfs { Vermes fc Kawkal 120081 ). Within 

the framework of the fossil field hypothesis this has been ex- 
Magnetic white dwarfs (MWDs) — white dwarfs with held n . n „ , . , rj , 

\ ' r; — ; -i| 3\ plained as a result of the average higher mass of the magnetic 

strengths ranging from about f kG (IJordan ct al. 2007) up ., . /T> TT , , ., , , ,. 
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6 , ITT — ; ,11 J J — fl progenitor systems Ap/Bp. However, detailed population 

to approximately f GG (Kawka ct al. 2007: K ulebi et al.l . , ,. , ', , , , , n , r » ;n 

I -f , : t* . — ~ , ; ^ synthesis studies show that the observed number ol Ap/Bp 

2009) — comprise more than about 10% of all white dwarfs. . . ~ . . . , ,, mm ,' , 

stars is insufficient to expfam the number of MWDs, and the 



There are two possibilities to account for the observed mag- 



contribution of an unseen pop ulation of weakly magnetized 



netic fields. According to the fossil field hypothesis, these . /T , . , , , , , \,, T . , — t~ 7^ n 

. 1 v ~ -11 1 A/B stars needs to be invoked (IW lckramasmgh c fc Ferrariol 

white dwarfs descend from Ap/Bp stars (lAngcl et al. 19811), UrUrh mi 1 -i -v^ ■ n . , ■ , n u , ra m 

„ fit- h- — 1= 11 iv> 2005). The second possibility is that high-field MWDs are 

and their magnetic helds are remnants of the previous cvolu- , , . ^ „ , , 

° produced m the aftermath of the merger of two mitiaify 

tion. One of the key properties of the population of MWDs is , . , .. , „ „, . , , 

fess massive white dwarfs, this scenario has been sug- 

their mass distribution, since it turns out that their average , , ... , . „ , . 

, , ' , , . , „ , , , gested as a viable mechanism for producing ultra-massive 

mass is considerably larger than that of the held popula- , ., , r 1 , N , , , . 

ITT — m _ , . _ white dwarfs (A7wd > I.IMq , and naturally explains 

tion (IKawka et al.H20071 ). Moreover, the reverse is also true. , n„, rT ^ • ,, ,i ■ n , , 

, : ~. — why MWDs are more massive than their held countcr- 

and the incidence of magnetism is larger m the population ™_. . „ . . , , r , , . 

parts, this scenario was hrst invoked for addressing the 

hot and massive population of white dwarfs in the ROSAT 

survey (|Marsh et al.l ll997T ). and recently to account for 
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|Kepler et all 120071 ) . Since the incidence of magnetism is 
higher for this population of ultr a-massive white dwarfs 
(|Wickramasinghe fc Ferrariol |2000| ). the merger scenario 
provides a natural expla nation for the properties of several 
MWDs, RE J 0317-853 (iBarstow et al.lll995l ; iFerrario et all 
Il997l: IVennes et all 12003 ) and SPSS J150746. 80+520958.0 
( Dobbie et al.ll201lh being some representative examples. 

There have been numerous suggestions that binary evo- 
lution could be respons i ble for mag netism in white d warfs . 
Specifically, iTout et all (|2008l ) and iNordhaus et all (|201ll ) 
discussed the possibility that turbulence during a common 
envelope phase could generate strong magnetic fields. These 
studies were aimed at explaining two observational facts. 
The first one is that no high-field MWDs have been detected 
in de tached binary systems with M dwarfs (|Liebert et all 
120051 ), while the second one is the higher incidence of mag- 
netism in white dwarfs in cataclysmic variables (approxi- 
mately 25%) with respect to tha t of sing l e star s (approxi- 
mately 10 %) — see, f o r inst ance, IWarner |l995). All these 
motivated ITout et all (|2008l ) to put forward the hypoth- 
esis that magnetism is a direct consequence of binarity, 
and specifically of the evolution during the co mmon enve- 
lope p hase of a close binary system. Specifically, I Tout et all 
(|2008l ) argued that isolated MWDs have merged, whereas 
those in magnetic cataclysmic variables have not. They 
therefore argued that the fields must have been generated 
before the merger, during the common envelope phase. Fol- 
lowing the merger there would be a period of normal red 
giant evolution in which the star would spin down until the 
magnetized core would emerge as a slowly rotating MWD. 
Only when the merger occurs just when the envelope is 
ejected would the whi te dwarf emerge rapid ly rotating. How- 
ever, a few years later iPotter fc Tout! (|20ld ) showed that the 
magnetic fields generated in this way are not durable if only 
simple diffusion of the field into the degenerate white dwarf 
matter is considered, because of the long associated time 
scales. 

There are other processes that might produce magnetic 
fields during binary evolution. One of such process is the ef- 
fect of shear during the intense accretion phase before igni- 
tion in an accreting white dwarf, a nd also dur ing the carbon 
simmering phase after the ignition jPirpll2008l) . By us i ng the 
physical reasoning of IPiro fe Bildstenl |20ol ). IPirol (|2008l ) 
proposed that poloidal and toroidal fields are produced 
due to a Taylor-Spruit dynamo during the differentially 
rotating phase. Their calculations show that for a white 
dwarf with Mwd = 1.37M , fi W D = 0.m K = 0.67 s -1 , 
M — 10 _7 MQyr~ 1 — where S7k is the Keplerian velocity 
— the shear layer produces a magnetic field with compo- 
nents B ^1O 8 G and B r w 10 5 G. 

The merger of two white dwarfs has become an ac- 
tive area of research during the last years. This is mostly 
due to the belief that double degenerate mergers might 
be successful progenitors of Type l a supernovae. Although 
this hypothesi s is relatively old (|lben fc Tutukovl 1 19841 ; 
IWebbink|[l984l ). this scenario has rec ently received a vig - 
orized interest — see, for instance, iBloom et all (|2012i ) 
and ISchaefer fc Pagnottal (|2012l) . Other reasons of in- 
terest in these models are the possib ilities that white 
dwarf mergers co uld form magnet ars (|King et all l200ll ; 
iLevan et all |2006|) and the dete c tion o f dust discs around 
white dwarfs (|Kilic fc Redfieldl 120071 ; | Garcfa-Berro et all 



20071; von Hippcl et al 



Brinkworth et al 
(iReid et~ 



2009) 



2001; 



120071 ; iGansicke et all |2007| . 120081 ; 

and dust discs around MWDs 



Dufour et all 120061 ; iFarihi et all 1201 ll ; 



iKawka fc Vennesfeoiil )" which might be indicative of a for- 
mer viscous disc phase following a merger process. 

In essence, modern hydrodynamical simulations of the 
merger of two white dwarfs show that the less massive sec- 
ondary is totally disrupted. Depending on the mass ratio 
of the binary components this process occurs on a dynam- 
ical time scale. As a result of the disruption of the sec- 
ondary, part of its material is accreted on the primary, 
whereas the rest of it forms a rapidly r otating Keplerian disc. 
This confirms the early sugg estions (|Tutukov fc Yungelsonl 
Il979l ; iNomoto fc ibenl I 1985) that a thick disc should be 
produced as an intermediary step following the coales- 
cence in which high accretion rates could be avoided and 
a SNIa explosion could take place. All these hypothe- 
ses have been recently confirmed b y the in-dep t h SPH 
simulations of iGuerrero et all (I2004T). lYoon et all (|2007l) . 
lLoren-Aguilar et al.l (|2009h . and fPan et all (|2012t ). In partic- 
ular, all these simulations show that a thermally supported 
shock heated region forms around the central object (the 
undisrupted primary) and, additionally, they also predict 
the formation of a massive thick disc supported by Keple- 
rian rotation. 

The accretion phase after the initial violent merger 
phase has also been subsequently investigated. Re- 
cent calc ulations of th e long-t erm evolution o f the 
merger dPiersanti et al. l2003bl lal; lYoon et all 120071 ; 
van Kerkwiik et al.ll2010h show that a supernova explosion 
is a likely outcome, although in the case of off-centre carbon 
ignition the result could also be an ac cretion induced 
collap se (AIC) to form a neutron star (|Nomoto fc Ibenl 
Il985h . However, there are multiple processes that influence 
the long-term evolution of the merger product. One of these 
is the int eraction of the white dwarf with the thin disc. 
Recently, ISaio fc Nomotol (|2004T ) studied this interaction 
th rough the v iscous bou ndary layer relying on the m odels 
of IPaczvnskH (|l99ll ) and IPopham fc Naravanl (|l99ll ). The 
current understanding is that an initial accretion stage 
nearly at the Eddington rate (about 10 -5 Mq yr _1 ) is 
unavoidable. More recently, another model has been pro- 
posed in which the thick disc evolves towards a spherical, 
thermally supported envelope due to the shear induced by 
the differential rotation within a sho rt timescale of 10 4 s 
(|Shen et all |2012| ; ISchwab et ail |2012| ). In both scenarios, 
the formation of ONe white dwarfs or, depending on the 
total mass of the system, of neutron stars via AIC, are the 
expected results. 

Finally, iGarci'a-Berro et all (|2012h proposed that the 
magnetic fields are generated in the corona above the merger 
product through an auj dynamo, since convection and dif- 
ferential rotation exist simultaneously. Using the results of 
the simulations of lLoren- Aguilar et all (|2009l ) they showed 
that very large magnetic fields can be generated in the hot 
corona resulting from a merger. Moreover, these magnetic 
fields do not diffuse into the white dwarf core or the sur- 
rounding disc, and remain confined to the outer layers, with 
strengths comparable to those observed in MWDs. Th i s find - 
ing has been recently corroborated by I Schwab et all (|2012l ) 
using the H0iland criterion. All in all, this means that mag- 
netic fields can be naturally produced in the merger of two 
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white dwarfs. However, one problem of this scenario is the 
excess of the remnant angular mo mentum, that is incompat - 
ible with observations of MWDs. iGarcfa-Berro et alj (|2012T ) 
suggested magneto-dipole radiation as a possible solution 
for this drawback of the model. However, it is likely that 
magneto-dipole radiation would not be efficient enough to 
slow down the remnants of the coalescence to the observed 
rotation periods of MWDs. Nevertheless, if the merger prod- 
uct has high magnetic dipole fields and the disc from the 
merger survives, it is possible that the long-term evolution 
of the system will involve angular momentum exchange be- 
tween the disc and the star, depending on the magnetic field 
strength of the central object. 

In this paper we analyse the long-term evolution of the 
disc and of the central object resulting from the coalescence 
of t wo white dwarfs, carrying over some e lements from mod - 
els (|Chatteriee et all l200d ; lAlparl l200ll ; lErtan et all 120091 ) 
constructe d for supernova fa llback discs around young neu- 
tron stars (jWang et al.l l 200(1 ). Our paper is organized as fol- 
lows. In Sect. 2 we explain the details of our model for the 
debris disc plus MWD s ystem, in Sect. 3 we u s e the results of 
the SPH simulations of lLoren-Aguilar et al.l (|2009l ) and the 
fallback model presented in Sect. 2 to compare our simula- 
tions with the observations of MWDs. Finally, in Sect. 4, we 
summarize our main findings, we discuss their significance 
and we draw our conclusions. 



2 THE DEBRIS DISC 

We consider a freely expanding Keplerian disc around the 
central object produced during the merger. This disc com- 
prises the mass of the disrupted low-mass companion that 
has not been incorporated to the central object, as little 
mass is ejected from the system and carries the excess an- 
gular momentum of the progenitor binary system. 



2.1 Evolution of the disc 

As mentioned, for debris discs the outer boundary of the 
disc is free. That is, it is not tidally truncated like discs in 
binary systems but, instead, it may expand. Such discs are 
similar to discs form ed by tidal disruption of stars by super- 
massive black holes ijCannizzo et al.ll 1990h or supernova fall- 
back discs around young neutron stars, e.g. a nomalous X- 
ray pulsars l|Chatteriee et alj|200d : lAlparll20o"ll ) . The model 
we em ploy is similar to the model proposed bv lErtan et al.l 
(|2009h which describes the evolution of the thin disc in terms 
of its initial mass and angular momentum. This description 
of the evolution of the disc presumes that the cause of an- 
gula r momentum transport in the disc is turbulent viscos- 
ity JShakura fe Sunvaev 19731) r ather than hydromagnetic 
winds (|Blandford fc Pavndl 19821 ) . 



The temporal evolution of the surface mass density, 
S, in the disc is described by a diffusion equation l|Pringlel 
[l98l]): 



<9E 
dt 



3d_ 

r dr 



.1/2 



dr 



vY>r 



1/2 



(1) 



where S is the surface mass density and v is the turbu- 
lent kinem atic viscosity . The rest of the thin disc structure 
equations (|Frank et al.l 120021 ) can be solved algebraically 



(|Cannizzo et al.l Il990l ) to obtain a viscosity in the form 
v = Cr p T, q where C, p and q are constant s deter mined by 
the opacity regime — see also lErtan et al.l (|2009l ) . For this 
viscosity Eq. CO is non - linear and has the self-similar solu- 
tions found by Pringld (|l974l ). In the first type of solution 
l)Cannizzo et al.l 1 1990h the angular momentum of the disc 
is constant (J c \ — 0), but the mass of the disc evolves as 
a power-law in time. Specifically, the mass flow rate at the 
inner boundary can be expressed as: 



M = M | 1 + — 

to 



M = (a-1)^ (2) 
to 



where to is the viscous time-scale at some characteristic 
radius, Mo the initial mass of the disc and Mo is the ini- 
tial accretion rate. The power-law index a can be written 
in term s of p and q, which a re determined by the opacity 
regime (|Cannizzo et al.l |l990i ) . For electron scattering the 
exponent is a = 19/16 and for the bound-free opacity regime 
is a — 5/4. Hence, the difference between the two regimes 
is not large, and similar results can be found by changing 
the initial mass of the disc slightly. For our numerical sim- 
ulations we use the bound-free regime as this would be the 
dominant opacity source for most of the lifetime of the sys- 
tem. 

The second type of solutions describes a disc in which 
the mass of the disc is constant (M<j = 0) while the an- 
gular mome ntum of the d isc increases by the torque acted 
by the star |Pringlelll99ll ). Such a disc would be relevant 
if the magnetosphere is rotating faster than the Keplerian 
velocity at the inner radi us of the disc. This occurs in the 
so-called propeller regime (|lllarionov fc Sunvaevlll975l ) if the 
mass is retained in the disc rather than being ejected from 
the system. 

The most important parameter governing the spin evo- 
lution of MWDs is the viscous time-scale to, which is deter- 
mined by the initia l mass and angular momentum of the disc 
(|Ertan et al.l l2009h . For our purposes, we assume a metal 
disc in the bound-free opacity regime with ass = 0.1, where 
ass is the Shakura-Sunyaev viscosity, for which the viscous 
time-scale for the white dwarf debris disc is 



to ~ 15 s 



3o 



10 18 cm 2 s" 1 



25/7 



Mo 



10- 1 MfT 







-3/7 



(3) 



where jo = Jo /Mo is the total specific angular momentum. 
Note that the viscous time-scale depends strongly on jo and 
also on the opacity. However, if the source of opacity is elec- 
tron scattering instead of bound-free transitions the viscous 
time-scale is changed at most by a factor of two. For a spe- 
cific angular momentum one order of magnitude larger the 
viscous time-scale is of the order of 1 day, which is compa- 
rable to the 3—14 days estimate of iMochkovitch fc Livid 
(|l99d ). 



2.2 Angular momentum evolution of the white 
dwarf 

Angular momentum exchange between a disc and a star has 
been extensively studied since the discovery of X-ray pul- 
sars, and these findings are nowadays applied to many as- 
troph ysical objects, like T Tauri stars and cataclysmic vari- 
ables (|Warnerfl995l) . The strong magnetic field of the white 
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dwarf disrupts the inner part of the disc which would oth- 
erwise extend to the surface of the star. The inner radius 
of the disc, R ln , is determined by the balance between the 
magnetic and the material stresses. This is of the order of 
the Alfven radius 



Ra 



= (- 



2/7 



\V2GMM J 



(4) 



which, in the case of spherical accretion, is determined 
by the balance between magnetic and ram pressures 
ijDavidson fc Ostrikerl 19731 ). Here G is the gravitational con- 
stant, /j,, is the magnetic moment of the white dwarf and M 
is the accretion rate reaching the the inner radius but not 
necessarily accreting. Thus, we use Ri n = £Ra where £ is a 
dimensionless factor used for converting the Alfven radius 
for spherical accretion to t he inner radiu s in disc accretion. 
In most studies — see e.g.. IWarnerl (|l995h — £ « 0.5 is usu- 
ally adopted. Furthermore, numerical studies have shown 
that this value is applicable to the case of a dipolar field 
(|Long et al.ll2005h . Nevertheless, when the accretion rate is 
so large that _Ri n is smaller than the radius of the star, we 
assume i?; n = Rwo- 

The way in which the star interacts with the disc de- 
pends on the interplay between its inner radius, the coro- 
tation radius, R co = (GM/D^-q) 1 ^ 3 where the disc is ro- 
tating with the same speed as the star, and the light cylin- 
der radius, _Rl = c/JIwd, c being the speed of light. To 
model it, it is customary to define the fastness parameter 
|Ghosh fe Lamblll979l ; [Lamblll988l . ll989l 1 as 



fi K (Rin) 



Rin 
Rco 



3/2 



(5) 



where Ok (-Rin) is the Keplerian angular velocity at the in- 
ner radius of the disc and f2wD is the angular velocity of 
the star. If the inner radius of the disc is smaller than the 
corotation radius (oj, < 1) the system is in the accretion 
regime. The inflowing matter at the inner edge of the disc 
will reach the surface of the star channelled by the magnetic 
field lines. If the inner disc radius is larger than the corota- 
tion radius (cj„ > 1) accretion cannot occur due to the fast 
rotation of the star. This cor r esponds to the prope l ler re gime 
jlllarionov fc Sunvaevl [l975l ; I Wang fc Robertsonl 1 19851 ). In 
this phase the disc continues its viscous evolution, though 
matter from the disc cannot be accreted on the star due to 
the centrifugal barrier. 

The torque appli ed by the disc to the white dwarf is 
jChosh fc Lamblll979l ): 



N d = n(bj.)N . 



(6) 



where No = (Rin) Rf n M and n ( cj«) i s call ed the di- 
mensi onless torque. Following lAlparl (|200lT ) and lEksi etHI 
(120051 1 we employ n = 1 — u>». Consider a ring of mass 
Am right outside the magnetospheric boundary 7?i n + AR 
where the flow is Keplerian. The angular momentum of 
the ring is Am(Ri n + AR) 2 flK(Rin + AR). Interaction with 
the magnetosphere will bring the matter to the angular 
velocity of the star at the magnetospheric boundary. The 
change in the angular momentum of the material is then 
Al = Am(R in + AR) 2 Q K (R in + AR) - AmR 2 n n t . Assuming 
the transition region is narrow (AR <C Ri n ) and a continu- 
ous mass flow we obtain N — MRf n Q,K(Rin)0- ~ LO *)- Hence 



this prescription is equivalent to assuming that the magne- 
tosphere and the accreting material are colliding "particles" 
and the interaction between them is inelastic. 

If the inner radius of the disc goes beyond the light 
cylinder radius, the disc can no longer torque the star via 
the magnetosphere, iVd = 0. In this so-called ejector stage 
the star will spin-down via the magnetic dipole radiation 
torque 



Amdr 



2(4 sin 2 gn|yp 
3c 3 



(-) 



where /3 is the inclination angle between the magnetic and 
rotation axis. 

We consider the effect of change in mass and radius 
of the star, which may be significant for those stars with 
masses close to the Chandrasekhar limit. In particular, the 
rotational velocity of the star depends on the moment of 
inertia of the star, which changes when a substantial amount 
of mass (of the order of about 0.1 M ) is accreted from 
the disc. In the case of accretion near the Chandrasekhar 
mass limit, the change in radius cannot be approximated by 
simple dependencies. Hence in our ca lculations we used the 
analytical mass-radius relationship of iNauenberd (| 19721) 



R* = 0.0126 R 



R 

R 



M 
M© 



-1/3 



i-(JL) 

\McnJ 



4/3' 



1/2 



(8) 



where 



Ch 



1.435 Mi 



■ 



(0) 



An important effect of the change in radius is that due to 
flux conservation the strength of the magnetic field changes, 
$ = BwdB|i D . This is also considered in our model. For 
the moment of inertia of white dwarfs we used 

0.34158 



h =3.2 x 10 50 gcm 2 



1.437 M 



© 



Mm 

M 

1.25" 



(10) 



which we obtained by fitting numerical results. 

Finally, we introduce a critical temperature below which 
turbulence in the disc can not be sustained. We find the 
effective temperature using 



(3GMM\ 



1/4 



Rir 



1/2 



1/4 



(11) 



(|Frank et al.ll2002l ) where n = N/My/GMRTn is the dimen- 
sionless torque per unit mass accretion rate acted on the 
star. 



3 RESULTS 

3.1 Discs from SPH simulations 

To compare our model with observations, we simulated the 
evolution of the merger remnants (namely, the white dwarf 
and the disc) using the results of the SPH simulations of 
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Table 1. Initial values o f our simulations, as obtained from 
lLoren-Aguilar et all j2009f ). 



Run 


M WD 


M 


SlwD 


Jo 


{M Q +M Q ) 


(M ) 


(Mq) 




(gems ) 


0.5+0.3 


0.62 


0.18 


0.134 


1.40 X 10 50 


0.8+0.4 


0.92 


0.28 


0.244 


2.65 x 10 50 


0.8+0.6 


1.10 


0.30 


0.401 


2.85 x 10 50 



lLoren-Aguilar et all (|2009h . We restricted ourselves to those 
simulations for which the total mass of the system does not 
exceed Mch- The initial conditions used in our simulations 
— that is, the masses of the central remnant and the disc, 
the angular velocity of the uniformly rotating star and the 
angular momentum of the disc — are shown in Table Q] 

For a given simulation, we calculated the spin evolution 
of the MWD for different polar field strengths. To determine 
the tot a l inte gration time we used the simple cooling law of 
iMestell (1 19651 ) and the otherwise typical observed value of 
T e ff«40 000 K. The exact value of of the initial core tem- 
perature depends on the mass of the merging components. 
However, its exact value has virtually no consequences in 
the age estimate, as white dwarfs cool very rapidly during 
the early phases of evolution. Consequently, we adopted an 
initial temperature of 10 8 K, a representative value of the 
SPH simulations. The cooling ages obtained in this way 
for the 0.5+0.3 M , 0.8+0.4 M and 0.8+0.6 M merg- 
ers — which result in MWDs of masses 0.8, 1.2, 1.4 Mq — 
were 13, 50 and 100 Myr respectively. It could be argued 
that for high core temperatures neu trino emission woul d 
shorten significantly the cooling times (|Althaus et al.ll2007h . 
However, given the simplicity of the approach adopted here 
we consider that these cooling times are accurate enough. 
Moreover, we found that when these cooling times are com- 
pared with_those_obtai^ calcula- 
tions (|Althaus et al. 2007; Sal aris et al.ll2O10h the differences 
are minor. 

In accordance with that explained in Sec. 12.21 the spin 
evolution of the star is determined by the fastness param- 
eter, to*. Initially, the accretion rate is very large and the 
magnetospheric radius is small. In this case, the fastness 
parameter is smaller than 1. Thus, the inner radius of the 
disc is equal to the stellar radius and the rapid accretion 
of matter from the disc increases the angular momentum 
of the star. Hence, initially the remnants of mergers are 
fast rotators, and they spin at almost one third of the crit- 
ical angular velocity of the star oj* « 0.33. This brief spin- 
up stage continues till u, > 1. At this point the magnetic 
torque from the disc starts to spin down the star. At this 
stage the star can reach critical angular velocities if the 
mass of the white dwarf is close to Mch- This is due to 
the steep dependence of radius on mass near Mch- Since 
spin-down is only effective for w* > 1, for very large accre- 
tion rates it might be possible that the star reaches critical 
rotation before the magnetospheric radius becomes larger 
than the stellar radius. This wa s the case for t h e simu - 
lations of IPiersanti et~al] (1200301). IPiersanti et all (|2003aT >. 
ISaio fc Nomotol (|2004l )7and lYoon et all (|2007l ). In our cal- 
culations this situation is avoided for the 0.5+0.3 Mq and 
the 0.8+0.4 Mq mergers. Hence, in the rest of the paper 



we only discuss these calculations. Additionally, it could be 
possible that the star might spin-up to a point in which the 
inner radius exceeds the light cylinder and, thus, accretion 
ceases. In this case, magnetodipole spin-down (a less effi- 
cient mechanism) takes over. However, in our calculations 
we do not find this to be the case. 

As the accretion rate decreases, the magnetospheric ra- 
dius, hence w*, increases rapidly. This results in a stronger 
torque on the star which eventually forces it to reach equi- 
librium at ~ 1. All our simulations evolve across this 
phase of torque transfer, and for all three cases there is 
a limiting magnetic field strength, Bij m , which determines 
whether a MWD will enter this strong spin-down phase, 
given the same cooling age. If the dipolar field strength, B p , 
is close to Bii m the final period of the star is very sensitive 
to the exact value of B p . In particular, near Pn m a 10% 
difference results in change by a factor of approximately 2 
in the final spin period. Finally, after this strong spin down 
phase, if the magnetic field strengths and cooling age are ad- 
equate, the system reaches quasi-equilibrium around w* < 1. 
This phase is called the tracking phase. The reason why a 
quasi-equilibrium phase is reached is that as mass transfer 
continues, the magnetospheric radius varies to match the 
corotation radius of the star. Hence, the star never reaches 
complete equilibrium and continues spinning down. 

The aforementioned evolutionary phases are common 
for the 0.5+0.3 Mq and the 0.8+0.4 Mq simulations, and 
the final period of the object is determined by B p , because 
in combination with the mass accretion rate, it determines 
the magnetospheric stopping radius, the fastness parameter, 
and thus the strength of the torque on the star for a given 
accretion rate. Since the accretion rate varies with time, the 
magnetic field strength of the star determines how fast the 
tracking phase is reached. In particular, if the dipolar field 
strength is large the tracking phase is reached earlier, be- 
cause the associated torques are larger. 

For the very early phases of evolution right after the 
merger, super-Eddington accretion rates are expected. This 
phase of the evolution is short (10 4 yr) with respect to 
the evolutionary time of the observed population of MWDs 
(10 7 — 10 s yr), but it determines the maximum angular ve- 
locities attained at the beginning of disc evolution, due to 
the initial spin-up phase. However, for a given white dwarf 
and disc angular momentum, the spin at which the tracking 
phase is reached is determined by B p . Actually, a simple re- 
lationship between B p and the final period of the star can be 
easily computed. Taking into account that during the track- 
ing phase the angular velocity of the star is almost equal to 
the Keplerian velocity at the inner radius of the disc, and 
the accretion is modelled by a power law, the period of the 
star can be written explicitly as: 

-10/14 , ^ ,6/7 



P, = 3.4 x 10' 



1.3 M 



R, 



4.3 x 10- 2 R 





18/7 



VHFG 
Jo 



2.5 x lO^gcm^- 1 



7.8 x 10 7 yr 



5/4 



25/28 



(12) 



which holds for B p > 10 MG. The results of our calculations 
are shown in Fig. [1] As can be seen, although the periods 
obtained in this way are more or less consistent with the 



© 0000 RAS, MNRAS 000, 000-000 



6 B. Kiilebi et al. 




Figure 1. Period versus magnetic field strength for MWDs with 
known periods, compared to our results. The solid line corre- 
sponds to the case in which a 1.2 Mq MWD is considered, while 
the dashed line corresponds to the case of a 0.8 Mm star. The 
known MWDs llKawka et al.1 12007| ; iBrinkworth et all 120071) are 
represented using black circles. The mass of the object deter- 
mines the size of the circle. The sizes of the circles for 1.3 Mq 
and 0.7 Mq stars are shown in the upper left corner of the fig- 
ure. The stars with unknown masses are represented using small 
black dots. The data points with arrows indicate lower limits for 
the respective periods. 



majority of periods of the populat ion of MWDs , which have 
periods of days to weeks — see iKawka et all (|2007n and 
references therein — the majority of the observed spin pe- 
riods seem to have a fairly fiat distribution between 1 and 
100 hr irrespective of the magnetic field. On the contrary, 
the models show a strong correlation with field strength in 
this range of periods. Before discussing this point we want 
to emphasize that not all magnetic white dwarfs have to be 
the result of a merger, and that our model is only aimed at 
reproducing the observed periods of massive white dwarfs. 
We elaborate on this point below. 



3.2 Low-mass discs 

Up to this point we have used an initial configuration di- 
rectly derived from the re sults of the SPH simulations of 
lLoren-Aguilar et al.1 (|2009h . Consequently, we have assumed 
that during the initial phases of the evolution the central 
white dwarf accretes at the Eddington rate, in accordance 
with previous works — see Sect. [[] However, an alternative 
model recently proposed bv lShen et ail (|2012h assumes that 
the disc disperses its angular momentum due to shear in the 
differe ntially rotating layers. The reasoning of IShen et al.1 
|2012l ) is based on the fact that radiative cooling is very inef- 
ficient. Hence, the material of the disc cannot relax to a thin 
disc configuration. This leads to entropy deposition within 
the disc due to viscous evoluti on. The detail e d mod elling of 
this phase was later done by ISchwab et al.1 (|2012h leading 
to viscous time scales of 10 4 s, within which the thick disc 
disperses and becomes thermally supported and spherically 



3 



cJj 




l°g(B,)/G 

Figure 2. Comparison of the calculations in which the results of 
SPH simulations are adopted (solid line), and the simulations in 
which a low-mass disc was adopted (black dots), for a 1.2 Mq 
MWD. The error bars represent the results of the simulations 
with different disc angular momenta. The upper error bars are 
the periods obtained when the disc angular momentum is half 
of that resulting from the SPH simulations, whereas the lower 
periods result from simulations in which the angular momentum 
is a factor of two larger than that of the SPH calculations. 

symmetric. Actually. IShen et al.1 (|2012l ) find: 



tvisc 



1 (r cy i\< 
cess \ h ) 



GM 



1/2 



3 X 10* s 



(13) 



where h is the thickness of the disc, and r cy i is its size. 

This model, however, has two drawbacks. The first one 
is that the model assumes that no accretion occurs what- 
soever. Consequently, the evolution of the disc is only gov- 
erned by the change in angular momentum — see Sect. 12.11 
This, in turn, forces the disc to evolve changing its angu- 
lar momentum, as opposed to the case where the mass is 
allowed to be accreted on the central object, in which case 
the disc can expand to conserve its angular momentum. The 
second drawback of this model is the assumption that radia- 
tive energy losses within the disc are negligible. A very crude 
approximation using an electron scattering opacity and the 
initial surface density of the disc leads to 



2 x 10 J s 



(14) 



where k is the opacity of the disc. This time-scale is or- 
ders of magnitude longer than t v i SC . However, both t v i sc and 
£rl depend strongly on r cy i. As a matter of fact it turns 
out that tRi/tvisc oc r cy i 1//2 . Hence, if the disc size increases 
during the evolution due to accretion on the central object, 
radiative cooling could become efficient and the disc would 
become thin. In our models, a thin disc reaches this size 
after 10 2 s if electron scattering is the dominant source of 
opacity. This is two orders of magnitude shorter than t v i BC . 
It is, nevertheless, important to realize that up to now we 
have adopted a power law for M. This approximation is safe 
for the later stages of the evolution, w hen most of the m ass 
of the disc could be lost if the model of IShen et al.1 (|2012l ) is 
valid. Hence, to estimate the final periods of the MWDs we 
also model these low-mass discs. In particular, we run a sec- 
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Figure 3. The evolution of the period of a 1.30 Mq MWD 
with a disc of mass 10~ 3 Mq for different initial magnetic field 
strengths. 



3.3 Some insight from RE J 0317 853 and 
WD 1658+441 

To compare our results with the properties of two known 
ultra-massive MWDs which are supposed to be the result 
of a merger — namely RE J 0317-853 and WD 1658+441 
— we followed the evolution of a remnant with their ap- 
proximate mass, 1.30 Mq. For this simulation we adopted 
the following initial conditions. The initial angular mo- 
mentum was Jo = 2.5 x 10 50 g cm s _1 , whereas the 
initial angular velocity was Qk/3. We remind the reader 
that, according to that discussed in Sect. 13.21 the value of 
the initial angular momentum plays a limited role, whilst 
the initial value of the angular velocity does not have 
a large impact in our results as long as we only con- 
sider the values of the periods obtained during the track- 
ing phase. The observed period of RE J 0317-853 is 721 s 
|Barstow et al.|[l995l ; iFerrario et al.|[l997l ') , whilst in the case 
of WD 1658+441 no conclusive eviden ce of rotation has 
been found so far jSchmidt et al.lll992T l. However, for this 
last object the instrumental timing li mit of 8.4 h can b e 
adopted as a lower limit to the period (|Shtol' et alj|l99^ ). 



170 MG and 660 MG 


Ferrario et all 1997; Burleieh et al.l 


1 19991: IVennes et al.ll2002 


: IFerrario & Wickramasinghell2005r). 



ond set of simulations in which we adopted the same total 
mass of the system, but we assumed a disc mass 10" 3 Mq 
— the mass of the disc when the Eddington accretion phase 
finishes and accretion driven by the evolution of angular 
momentum takes over — while the rest of the material was 
incorporated into the remnant white dwarf. We performed 
this test for the 0.8+0.4 Mq merger — which in this case 
results in a MWD of about 1.2 Mq. Our results are shown 
in Fig. [2] As can be seen, the final periods of these simula- 
tions are almost identical to those previously presented. The 
reason for this is that the final periods reached at the end of 
the simulation for dipolar field strengths below Bii m are sen- 
sitive to the starting angular momentum of the star, since at 
this phase the star carries most of the angular momentum 
from the merger process. However, for stars with dipolar 
field strengths above Bn m , the final periods are determined 
by the angular momentum of the disc, see Eq. (|12[) . This 
equation is valid for cooling ages longer than about 10 7 yr 
and Bp > Bn m « 10 MG, which apply for most high-field 
MWDs with known periods. 

The angular momentum of the disc depends on the mass 
of the primary and the mass difference between the two com- 
ponents. As can be see n in Table [p the angular m omenta of 
the mergers studied bv lLoren-Aguilar et al.l (|2009l ') differ by 
about a factor of two. To simulate arbitrary mass ratios, we 
made trial runs with different disc angular momenta within 
a factor of 2 of those obtained in the SPH simulations — 
see Fig. [5] Our calculations show that, except for the case 
in which the magnetic field strengths are close to Bn m - 
that is, when the evolution is stopped just before the quasi- 
cquilibrium phase is reached — the results are nearly identi- 
cal to those in which the angular momentum resulting from 
the SPH simulations is adopted. In fact, although for the 
strong torque transfer phase the periods can differ by an 
order of magnitude, for the tracking phase the maximum 
difference in the computed periods is approximately 20%. 



the magneti c field strength of WD 1658+441 is estimated to 
be 2.3 MG (|Schmidt et al.lll992r >. For these specific simu- 
lations we stopped the simulations at the cooling age for 
which T c ff ~ 40 000 K, a value consistent with the obser- 
vational determination for RE J 0317—853 jBarstow et al.l 
1 19951 ; IFerrario et al.|[l997l ; IVennes et al.ll2003T l and 30 000 K, 
the effective temperature of WD 1658+441. However, we 
stress that if T c g = 40 000 K is taken as an upper limit, the 
final ages for both objects are lower limits {t coo \ « 10 s yr) 
and consequently the computed periods are also lower lim- 
its. The evolution of the period as a function of time is 
shown in Fig. [3] Clearly, Fig. [3] reveals that our simula- 
tions cannot reproduce the observed properties of these ob- 
jects, and that another braking mechanism must operate. 
In particular, for RE J 0317—853 our model predicts a pe- 
riod which is much longer than that observationally deter- 
mined, while for WD 1658+441, which has a smaller mag- 
netic field, we obtain periods of the order of 10 2 s, much 
shorter than expected. More generally, for dipole magnetic 
field strengths larger than Bii m ~ 10 MG — namely, for the 
case of RE J 0317—853 — the periods are longer than about 
10 4 s, while for field strengths smaller than this value — 
which is the case of WD 1658+441 — the spin-down mech- 
anisms implemented in our model are not adequate either. 
Hence, our simulations show that neither viscous dissipa- 
tion nor magneto-dipole spin-down are able to explain the 
observed periods of these objects, confirming our previous 
finding that other alternatives must be sought. 



3.4 Models with turbulence turn-off 

An important question to be analysed is whether turbulence 
in the accretion disc can be sustained efficiently during the 
entire lifetime of MWDs. This question is relevant when the 
lifetimes of neutron stars with fallback discs and MWDs 
are compared. In the former case the evolutionary ages are 
typically 10 4 — 10 s yr, whilst the lifetimes of MWDs are three 
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Figure 4. The evolution of the period for a 1.30 Mq MWD. 
The simulations are identical to those of Fig. [3] except for the 
inclusion of turbulence turn-off at T p = 1 000 K. The dashed lines 
correspond to the phase during which the torque transfer from 
the disc becomes ineffective and magneto-dipole spin down takes 



to four orders of magnitude longer. Thus, this question is an 
important factor to be considered. 

Whether the disc is turbulent or not is mainly deter- 
mined by the degree of ionisation of the disc material. That 
is, the existence of turbulence depends on the tempera- 
ture, which, in turn, is determined by the accretion rate 
— see Eq. (JTTJ) . The ionisation sta bility of accretion d iscs 
has been investigated thoroughly bv lMenou et~aH i|200ll ). In 
their work they computed the thermal equilibrium (the so- 
called S-curves) of thin accretion discs with different chem- 
ical compositions, and showed that for thin discs made of 
C, thermal equilibrium is reached near 5 000K, and that 
for all compositions recombination always occurs below ap- 
proximately 1 000 K. For the conditions found during the 
final stages of our simulations these temperatures are rel- 
atively high. However, it is known that these discs can be 
partially ionised at even lower temperatures, due to sev- 
eral other processes, like cosmic ray heating, coll isions and 
charged-particle mobility. Ilnutsuka fc San 

3 d200l l carefully 

studied all these effects and reached the conclusion that at 
temperatures below about 1 000 K, magneto-rotational in- 
stabili ty (MRI) co u ld be still effective. Building upon their 
work, lErtan et all {{200$ ) modelled the evolution of debris 
discs assuming that the differentially-rotating layers with 
temperatures smaller than a critical temperature, T p , do not 
change the angular momentum of the star. Moreover, they 
showed that when T p ~ 100 K is adopted the observations 

are consistently reproduced. 

Contrary to the modelling of lErtan et al.l (|2009h . in 
which turbulence is switched-off layer by layer, we applied 
a simpler procedure to test for how long the disk would be 
turbulent. For each time step, we calculated the tempera- 
ture at r = 10i?in. If this temperature was smaller than 
T p = 1000K, we assumed that the disc was not turbulent 
and accretion was switched-off. Thus, the MWD immedi- 
ately enters in the so-called ejector phase and spins down 
only due to magneto-dipole radiation. To compare these cal- 
culations to the case of RE J 0317-853 and WD 1658+441 



we applied this procedure to the model of 1.3 Mq described 
in Sec. 13.21 We found that although the MRI is long-lived 
( about 10 s — 10 7 yr), the disc becomes inactive at the cur- 
rent ages of these objects (about 10 s yr). Moreover, since 
the magneto-dipole torque is not as strong as the magnetic 
torque from the disc, the final period is determined by the 
time at which MRI is turned off, which depends on the ac- 
cretion rate, and can be determined from Eq. (Ill[) : 



P, = 1.3xl0 3 sf M * 



1.3 Mq 



AL 



-10/28 



1/28 



2.2 x 10" 11 Mgyr- 1 



lit, 



3/7 



10 s MG J V 4.3 x 10- 3 Rp 



© 



-9/7 



1000K 



(15) 



As can be seen in Fig.0] the accretion rate at the time of 
MRI deactivation depends also in the strength of the mag- 
netic field. If the torque exerted by the disc is strong enough, 
the disc becomes inactive after reaching the tracking phase. 



For T D 



1 000 K, for white dwarfs of masses between 0. 



and 1.3 Mq this occurs for magnetic fields around 90 MG. 
In this case, the final period of the star is determined by B p 
and the critical temperature of MRI turn-off. Using Eq. (JTTJ) 
once more, the value of M at the point at which T p is reached 
can be determined, and employing this value, the period at 
the time of MRI turn-off can be computed: 



P* = 1.3 x 10 3 s 



11, 



M* 



1.3 M 



-5/13 



4.3 x 10- 3 R, 



■0 



18/13 



10 8 G 



6/13 



1000K 



-12/13 



(16) 



This expression shows that for MWDs with field 
strengths between 50 MG and 1 GG, and T p ranging from 
300 K to 1 000 K, the final periods reached are around 
10 3 - 10 4 s (see also Fig.Hji. 

In Fig. [S] we compare the results of these calculations 
with the observational data. This is done for three white 
dwarfs masses, 0.8 Mq, 1.0 Mq and 1.3 Mq. As can be 
seen in this figure, now the computed periods are consider- 
ably shorter, and match better the periods of MWDs with 
known masses. However, in order to perform a meaningful 
comparison with the observations only those MWDs which 
are likely to be the result of a merger need to be consid- 
ered. Thus, we removed all MWDs with masses smaller than 
0.8 Mq. Of these MWDs six have known periods. Three of 
them are well known slow rotators, namely Grw +70°8247, 
G 240-72 and WD 1658+441, for which only lower limits to 
their periods have been determined. The remaining three are 
PG 1015+014, PG 1031+234 and RE J 0317-853, which all 
have magnetic fields above 50 MG, and periods of the order 
of an hour or less. We simulated the secular evolution of the 
spin of these three MWDs, using the observationally deter- 
mined values of their masses and magnetic field strengths — 
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Table 2. Comparison of the simulated and observational data, when a an MRI turn-off temperature T p = 1000K is adopted. 



WD 


Other names 


(M ) 


B 

(MG) 


Bp 
(MG) 


^obs 
(min) 


1 sim 

(min) 


Refs. 


0325-857 
1015+014 
1031+234 


RE J 0317-853 
PG 1015+014 
PG 1031+234 


1.35 
1.15 
0.93 


170-660 
50-90 
200-1000 


200 
90 
500 


12.08 
98 84 +0 - 14 
211.8 + 3.0 


28 
36 
205 


1,2,3 
4,5,6 
4,6,7 



m lBarstow et alj <1995ft. (2 ) IFerrario et al.l ||l997l). (3)lKiilebi et al.l bOlOh, (4HLiebert et al.l [120031) , 
(5) lEuchner et alj J2006T) . (6) iBrinkworth et alj j2007T ), (7) ISchmidt et alj Jl986h 
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Figure 5. Same as Fig.[T]for the simulations in which turbulence 
is switched-off at T p = 1 000 K. The solid line represents the spin 
evolution of a MWD of 1.3 Mq, the dashed line that of a MWD 
of mass 1.0 Mq, and dot-dashed one corresponds to the case in 
a mass of 0.8 Mgis adopted. 



see Tabled In addition to the properties of the disc, the most 
important input in our simulations is the dipolar magnetic 
field strength. Sp ecifically, it has been shown that probably 
RE J 0317-853 (iBarstow et alj [l995l; IFerrario et all Il997l ; 



IVennes et alj|2003l). PG 1031+234 jSchmidt et al.lll986l ) and 
PG 1015+014 (|Euchner et alj|2006j ) have dipolar field com- 
ponents considerably smaller than the maximum observed 
field strength. Thus, for them we adopt the values listed 
in the fifth column of Table [21 which in general are some- 
what smaller than the observed surface field strengths. The 
results of our simulations are shown in Table[2] We note 
that the simulated periods for these three MWDs (listed 
in the seventh column of this table) compare favourably 
with the observed periods (sixth column). In particular, 
for RE J 0317—853 the computed period is a factor of two 
longer, while for PG 1015+014 is a factor of three shorter, 
and for PG 1031+234 we find a nice agreement. All this in- 
dicates that these objects were most likely originated in a 
merger, and subsequently lost their angular momentum due 
to the magnetospheric interaction with a remnant disc. 

3.5 Winds 

It is currently understood that the initial rapid accretion 
phase can not be approximated by cold accretion, and that 
a hot re gion around the cen t ral M WD would be formed — 
see, e.g.. Ivan Kerkwiik et"aH (|2010h . In order to qualitatively 




Figure 6. The evolution of the period of a 1.1 Mq MWD due to 
torques from the wind with turbulence turn-off at T p = 1000K. 
The initial conditions were those obtained in the SPH simulation 
of the merger of two 0.8+0.6 Mgwhite dwarfs (see Table[TJl. 



understand the implications of this for the angular momen- 
tum evolution of a MWD, we built a model in which the 
spin down is due to the outflow of the material rather than 
to accretion onto the star. In this prescription, mass transfer 
between the disc and the star occurs, but the material is not 
bound to the system, and it is ejected. 

The induced torque s on the star can be then approxi- 
mated by wind torques (|Weber fc Davis!ll967l ): 



Ra 



(17) 



where M w is the wind mass loss rate, and 77 is a parameter 
which depends on the geometry of the mass loss, which for 
a spherically symmetric wind is 2/3. In this expression the 
Alfven radius Ra is computed as the distance at which the 
Alfven speed equals the wind velocity, rather than the in- 
coming velocity of the accreting material, as in Eq. (Q. Note 
that Ra depends on the geometry of the magnetic field, and 
simulat ions are necessary for precis e values — see, for in- 
stance, [MatTfePudriti J200!, [2008|). However, for the sake 
of simplicity, we assume that the mass-loss rate follows the 
same power law given by Eq. ([2]) and that Ra is that corre- 
sponding to disc accretion. 

We find that during the early phases of the evolution the 
spin rates are not strongly affected, since the mass-loss rates 
are very high, and Ra is small. However, as the mass-loss 
rate decreases Ra becomes larger, and a strong torque on 
the object is exerted. In fact, the most apparent qualitative 
difference between the wind and the disc spin down scenarios 
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is that in the first one the torque spins down the MWD most 
effectively at relatively early phases of its evolution, before 
reaching a quasi-equilibrium, whereas in the second one the 
wind spin down is ineffective during the initial phases of 
the evolution, and as the mass-loss rate drops, the torque 
rapidly becomes stronger. Moreover, in this case the period 
evolution can be calculated analytically, once it is taken into 
account that the only time-dependent quantity now is M: 



(18) 



In (JL) oc J±M* /7 Rl 4/7 M* /7 t 3 /28 B s/7 t ri/2S 



where Po is the initial period of the star. This expression 
indicates that the final periods of MWDs are very sensitive 
to the values of the stellar parameters, and that depending 
on their precise value large spin downs might be obtained. 

We are interested in estimating the time needed to reach 
periods similar to that of RE J 0317—853, i.e. P* « 10 4 s. 
Assuming that the initial period is about 10 s, the corre- 
sponding time-scale is: 
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For strong dipolar field strengths periods of the order of 
10 3 seconds are reached in time-scales much smaller than the 
cooling age of RE J 0317—853 and similar objects. Addition- 
ally, our simulations with turbulence turn-off show that for 
MWDs with field strengths larger than 200 MG the wind 
spin down is effective long enough to induce strong spin 
down — see Fig. [6] We also note that in this case the re- 
sulting periods are similar to those of Grw +70 °8247and 
GD 229, which > 100 yr — see lKawka et all (|2007f ) and ref- 
erences therein. On the contrary, for field strengths smaller 
than about 50 MG, the wind torques are ineffective, yielding 
periods smaller than a few minutes (P» < 3 x 10 2 s) which 
are unrealistic for the observed periods of the MWDs. 



4 CONCLUSIONS 

We have studied the long-term evolution of magnetic white 
dwarfs surrounded by debris discs, which are thought to be 
the final products of the merger of two degenerate cores. To 
do so we built a simple model to follow the coupling between 
the central object and the surrounding disc. This model 
is similar to those used to study supernova fallback discs 
around young neutron stars, and presumes that the evolu- 
tion of the disc is dominated by angular momentum trans- 
port by turbulent viscosity. Contrary to previous studies of 
the evolution of discs around white dwarfs, our model con- 
siders a freely expanding Keplerian disc. Using this model we 
found that for these systems there is a relation between the 
mass, the rotation period, and the strength of the magnetic 
field. 

In a first set of simulations we studied the evolution of 
the merger remnants using the results of detailed SPH sim- 
ulations and assuming that the entire disc was turbulent. 



For these models we found that the final spin of the white 
dwarf is given by a simple relation — see Eq. I)12fl • However, 
the periods obtained in this way are somewhat larger than 
the observed ones. Consequently, we studied other alterna- 
tives. In particular, in a second set of simulations we com- 
puted a set models in which we assumed — in accordance 
with the studies of lYoon et al.1 (120071 ); I van Kerkwiik et al.l 
l|2010h ; IShen et al.1 (|2012h — that the cold accretion as- 
sumption would not be applicable to the post mer g er sys - 
tems. Following the suggestion of iNordhaus et all (|201ll ). 
we simulated the case in which the mass transferred from 
the dies would end up as an outflow, and the spin down 
is governed by the magnetic interaction with it. We found 
that, nevertheless, the periods obtained in this second set 
of models are again too large. We also varied by a factor 
of two the initial angular momentum of the disc, and we 
also found that the resulting periods do not differ much 
from those obtained using the results of the SPH simula- 
tions. All this led us to analyse whether turbulence in the 
accretion disc can be sustained efficiently during the entire 
lifetime of MWDs. This depends mainly on the degree of 
ionisation of the disc material. According to the most re- 
cent studies we assumed that turbulence ceases when the 
temperature of the disc is smaller than a certain thresh- 
old, f or which we adopted T p = 1 000 K (|lnutsuka fc Sane] 
l2005h . When this model is adopted our calculations repro- 
duce the observed properties of three MWDs with known 
masses and periods in the solar neighbourhood which are 
massive enough — namely, with masses larger than 0.8 Mq 
— to be considered the result of white dwarf binary mergers 
(RE J 0317-853, PG 1015+014, and PC 1031+234). These 
findings support the hypothesis that these white dwarfs are 
the remnants of double-degenerate mergers that have gone 
through magnetopsheric interaction with a remnant disc. 
More specifically, it has been known for a long time that 
RE J 0317-853, the fastest rotator within the MWD pop- 
ulation, spins slower than what it should be expected for a 
double white dwarf merger. Our model reproduces the rota- 
tion period of this object, which was previously unaccounted 
for. However, we also find that our model cannot explain the 
properties of the rest of the population of massive MWDs, 
whi ch is probably formed by multipl e evolutionary chan- 
nels l^narig^^Wickrarnasinghe] (|2005l ). In particular, there 
are three massive MWDs (Grw +70°8247, G 240-72 and 
WD 1658+441) which are slow rotators, and probably have 
experienced spin-down during the red giant phase either as 
a single star of after a core m erger durin g the c ommon en- 
velope phase, as suggested bv lTout et all (|2008h . 

iFerrario fc Wickramasinghel (|2005l ) conducted very 
thorough observational investigations of the properties of 
the population of MWDs and resolved that within it there 
are three different magnetic field and period intervals, which 
point out towards different evolutionary origins. Specifically, 
they proposed the existence of three distinct groups of stars. 
The first group corresponds to a population of magnetized 
slow rotators — with typical periods ranging from P « 50 yr 
to 100 yr — which most likely have evolved from isolated 
progenitors. To the best of our knowledge only one qual- 
itative study of the evolution of rotating MWDs originat- 
ing from single st ars has been performed so far. Specifically, 
ISuiis et al.l (|2008f ) estimated the effect of magnetism on the 
core rotation and concluded that the presence of magnetic 
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fields leads to slower rotating cores. However, their results 
do not fully account for the long rotation periods of some 
MWDs, namely those belonging to the tail of the distribu- 
tion (with periods longer than about 100 yr). For these stars 
they proposed that spin-down by stellar winds might be an 
effective alternative. However, a binary origin fo r this group 
of sta rs cannot be discarded, as the results of I Tout et al.l 
|2008l ) show that a merger during the common envelope 
phase can result as well in a slowly rotating MWD. A sec- 
ond group o stars is that of strongly magnetized fast rota- 
tors (P ~ 700 s), which probably carry the remnant angu- 
lar momentum from the double-degenerate merger. Finally, 
there exists a third group of rotators with periods ranging 
from hours to days, which may have a mixed origin. Addi- 
tionally, they found that the majority of MWDs have peri- 
ods between a few hours to a few weeks, and also that the 
known slow rotators are not all necessarily highly magnetic 
(e.g. WD 1 658+441). However, althoug h the o verall picture 
drawn by iFerrario & Wickramasinghel (|2005h may not be 
significantly altered by our calculations, our results indicate 
that the distinction between all three groups may not be 
so clear — see Fig. [5] Moreover, our results show that slow 
rotation does not provide sufficient evidence to discard a bi- 
nary evolutionary origin of individual MWDs, since also the 
dipole field strength, the mass and the age of the specific 
star need to be considered. 

Additionally, our scenario predicts the existence of a 
population of MWDs surrounded by dust discs. Observa- 
tional searches of dust discs around two massive MWDs 
(RE J 0317-853 and WD 1658+441) have been performed 
using the Sp itzer telescope, but no discs have been de- 
tected so far (|Hansen et alj|2006h . However, the lack of de- 
tection of debris discs around high-field MWDs (those with 
B > 2MG) might be due to the diffuse structure resulting 
from their large magnetospheric stopping radii — a con- 
sequence of the large magnetic pressure of the remnant of 
the merger — which, in turn, produces an extended disc. 
More generally, our models imply that the detectability of 
a remnant disc depends not only on the mass of the MWD, 
but also on its magnetic field strength and on its period. 
Clearly, deeper observations are needed to confirm the exis- 
tence of these diffuse discs around MWDs. In this regard it 
is important to mention the recent discovery of an unveiled 
population of MW Ds surrounded by dust discs. NLTT 10480 
(|Kawka fe VennesH201ll l and G77-5018 |Farihi et al.ll201lh 
are examples of MWDs with significant infrared excesses, 
which could be attributed to the existence of debris disc 
around them. 

Finally, we would like to mention that one important 
hindrance of our modelling is the difference between the ge- 
ometry of the magnetic field assumed in our simulations and 
that of the observed MWD population. Multi-ph ase spectro- 
polarimetric analysis dEuchner et al.|p002). I2005T) and single 
phase spectral studies ( Kiilebi et al.l 12003 ) of a large num- 
ber of MWDs have showed that their field geometry is more 
complicated than a simple dipole, as it is assumed in our 
simplified model. This is an important issue as numerical 
simulations show that the geometry of the magnetic field 
is important for determinin g the magnetospheric s t oppin g 
radius of an accreting star (|Long et al1l2007l . l200St l2012h . 
For the case of RE J 0317-853 and PG 1031+234 — two 
of the specific white dwarfs analysed here — we addressed 



this issue by assuming dipolar field strengths smaller than 
the total measured value of the magnetic field, a simple but 
effective procedure. Thus, future studies in this direction 
should adopt more realistic magnetic configurations. 
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